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Abstract 

The present methods for obtaining the optimal Lewenestein- Sanpera decomposition 
of a mixed state arc difficult to handle analytically. We provide a simple analytical 
expression for the optimal Lewenstein-Sanpera decomposition by using semidefinite pro- 
gramming. Specially, we obtain the optimal Lewenstein-Sanpera decomposition for some 
examples such as: Bell decomposable state, Iso-concurrence state, generic two qubit state 
in Wootters's basis, 2 3 Bell decomposable state, d d Werner and isotropic states, 
a one parameter 3 (g) 3 state and finally multi partite isotropic state. Keywords: Op- 
timal Lewenstein-Sanpera decomposition, Semi-definite programming, Bell 
decomposable states, Werner and isotropic states. 

PACs Index: 03.65.Ud 
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Entanglement is one of the most striking features of quantum mechanics [1, 2]. In the case of 
pure states it is easy to check whether a given state is, or is not entangled. For mixed states, 
however, the statistical properties of the mixture can hide the quantum correlations embodied 
in the system, making thus the distinction between separable and entangled states enormously 
difficult. 

In the pioneering paper [3], a very interesting description of entanglement was achieved 
by defining the best separable approximation (BSA) of a mixed state. In the case of 2- 
qubit system, it consists of a decomposition of the state into a linear combination of mixed 
separable part and a pure entangled one. In this way, the whole non-separability properties 
are concentrated in the pure part. 

In the Ref. [3], the numerical method for finding the BSA has been reported. Some 
analytical results are also obtained for special states of two qubit states [4] . Further in [5] the 
BSA of two qubit state has been obtained algebraically. They have also shown that in some 
cases the weight of the entangled part in the decomposition is equal to the concurrence of the 
state. An attempt to generalize the results of Ref [3] is made in [6]. 

There is another method we can use which achieves exactly the same effect, called semidef- 
inite programming(SDP). Over the past years, SDP in particular, have come to be recognized 
as valuable numerical tools for control system analysis and design. In SDP one minimizes a 
linear function subject to the constraint that an affine combination of symmetric matrices is 
positive semidefinite. SDP, has been studied (under various names) as far back as the 1940s. 
Since 1990 many applications have been discovered in areas such as estimation, signal pro- 
cessing and it is currently considered to be the hottest area in optimization. Although SDP is 
designed to be applied in numerical methods it can be used for analytical computations. All 
of the above mentioned apphcations indicate, that the method of SDP is very useful. Some 
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authors try to use the SDP to construct an exphcit entanglement witness and entanglement 
distillation[7, 8, 9]. 

In this paper we use the SDP method in order to obtain the optimal Lewenstein-Sanpera 
decomposition (LSD) of a mixed state. Then we show how to perform the optimal LSD for 
well known different examples via SDP method. 

The paper is organized as follows: 
In section-2 we define SDP. In section -3 we give a brief review of Optimal LSD. In section 
-4, using the SDP method we obtain the optimal LSD of some mixed state density matrices 
such as: Bell decomposable state, Iso-concurrence state, generic two qubit state in Wootters's 
basis, 2 (8) 3 Bell decomposable state, d<Sid Werner and isotropic states, a one parameter 3 (8) 3 
state and finally multi partite isotropic state . The paper is ended with a brief conclusion. 

2 Semi-definite programming 

A SDP is a particular type of convex optimization problem [10]. A SDP problem requires 
minimizing a linear function subject to a linear matrix inequality (LMI) constraint [11]: 

minimize V = c^x 

(2-1) 

subject to F{x) > 0, 

where c is a given vector, = {xi, ...,Xn), and F{x) = Fo + J2i XiF^, for some fixed hermitian 
matrices Fj. The inequality sign in F{x) > means that F{x) is positive semidefinite. 

This problem is called the primal problem. Vectors x whose components are the variables 
of the problem and satisfy the constraint F{x) > are called primal feasible points, and if 
they satisfy F{x) > they are called strictly feasible points. The minimal objective value c^x 
is by convention denoted by V* and is called the primal optimal value. 

Due to the convexity of set of feasible points, SDP has a nice duality structure, with, the 
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associated dual program being: 



maximize 



■Tr[FoZ] 



Z >0 



(2-2) 



Tr[FiZ] = a. 



Here the variable is the real symmetric (or Hermitean) matrix Z, and the data c, Fj are 
the same as in the primal problem. Correspondingly, matrices Z satisfying the constraints 
are called dual feasible (or strictly dual feasible if Z > 0). The maximal objective value of 
—TtFqZ, i.e., the dual optimal value, is denoted by d*. 

The objective value of a primal(dual) feasible point is an upper (lower) bound on V*{d*. 
The main reason why one is interested in the dual problem is that one can prove that d* <V*, 
and under relatively mild assumptions, we can have V* — d*. If the equality holds, one can 
prove the following optimality condition on x: 

A primal feasible x and a dual feasible Z are optimal which is denoted by x and Z if and 



This latter condition is called the complementary slackness condition. 

In one way or another, numerical methods for solving SDP problems always exploit the 
inequahty d < d* < V* <V, where d and V are the objective values for any dual feasible 
point and primal feasible point, respectively. The difference 



is called the duality gap. If the equality d* — V* holds, i.e., the optimal duality gap is zero, 
then we say that strong duality holds. 



only if 



F{x)Z = ZF{x) = 0. 



(2-3) 



'p*-d*^ c^x + Tr[FoZ] = Tr[F{x)Z] > 



(2-4) 
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3 Optimal Lewenstein-Sanpera decomposition 

According to pioneering work of Lewenstein and Sanpera [3], any bipartite density matrix p 
has a decomposition of the form 



where p'^ is a separable density matrix, pe is a entangled state, and the parameter A e [0, 1]. 
According to the theorem (2) of reference [3], any 2— qubit density matrix p can be written as 



where the entangled part is pure state, but in general pe can be a mixed or a pure state where 
the whole entanglement of p is concentrated in (1 — X)pe- 

Rare exceptions aside, the LSD of a given (non separable) p is not unique, there is usually 
a continuum of LSD to choose from. The decomposition with the largest weight A of the 
separable part is the optimal LSD with respect to the chosen separable set, which is proved to 
be uniquely determined. According to its definition, the separable part of this decomposition 
is called the best separable approximation (BSA) of p, and its weight A the separability. 

4 Optimal Lewenstein-Sanpera decomposition with semi- 
definite programming 

In this section using the SDP method we obtain the optimal LSD of a mixed state for some 
well known different mixed states. 



p = (1 - X)pe + Xp, 



(3-1) 



P = Ap: + (1-A)|^)(^| 



(3-2) 
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4.1 Optimal Lewenstein-Sanpera decomposition for Bell-decomposable 
state 

A Bell decomposable (BD) state is defined by: 

4 4 

i=l 

where j^'i) is Bell state, given by: 



i=l 



V2 



(ITT) + IU)), 



r) = ;^(ITT)-|U)), 



1^4) = ^ 



^(in) + iiT)), 



V2 



(in) -lit))- 



(4-1) 



(4-2) 
(4-3) 
(4-4) 

(4-5) 



In terms of Pauli's matrices, p can be written as. 



1 ^ 

p= -(7{8)7-K^ij(7i(8)(7j), 

i=l 



where [12] 



^1 = P1-P2+P3-P4, 

^3 = Pl+P2-P3-P4- 



Prom the positivity of p we get 



(4-6) 



(4-7) 



l+ti-t2 + t3> 0, 

1 - tl + t2 + ts > 0, 

l+ti+t2-t3> 0, 

l-ti-t2-t3> 0. 



(4-8) 
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These equations form a tetrahedral with its vertices located at (1, —1, 1), (—1, 1, 1), (1, 1, —1), 
(—1, —1, —1) [12] . In fact these vertices denote the Bell states given in Eqs. (4-2) to (4-5), 
respectively. 

On the other hand p given in Eq. (4-6) is separable if and only if ti satisfy Eq. (4-8) 
together with the following equation 

l+ti+t2 + t3> 0, 
l-tl-t2 + t3> 0, 

(4-9) 

l + ti-t2-ts> 0, 
l-ti+t2-t3> 0. 

Inequalities (4-8) and (4-9) form an octahedral with its vertices located at Of = (±1, 0, 0), 
Ot — (0, ±1,0) and Of = (0,0, ±1). So, tetrahedral is divided into five regions. Central 
regions, defined by octahedral, are separable states {pk < |). There are also four smaller 
equivalent tetrahedral corresponding to entangled states (pjt > | for only one of A; = 1, ...,4), 
where Pk — \ denote to boundary between separable and entangled region. Each tetrahedral 
takes one Bell state as one of its vertices. 

Now in order to obtain optimal LSD of entangled BD state given in (4-1), with pi > |, we 
first choose an arbitrary separable state 

Ps = jlPi \<t>i) {4>i\. < < 1, j^Pi = 1, p'l < i jZPi = 1' (4-10) 

i=l i=l ^ i=l 

in the separable region. Then using strict SDP optimization prescription of section (2), we try 
to optimize Tr(Ap'^) with respect to p — Ap'^ > 0, where the feasible solution corresponds to 

^max = min{^,-^,-^,-^}. (4-11) 
Pi P2 Ps Pi 



Now, using the inequalities 



Amax<^, /ori = 2,3,4 (4-12) 

Pi 
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and summing over the indices i=2,3 and 4, we obtain 

{l-p[)Amax<{l-pi), (4-13) 

since, we have 

The only possible choice of A^ax which is consistent with the positivity of p — A^^xP's is 

Amax = rr^v (4-15) 

(1 - Pi) 

This choice of Aj^ax given in (4-15) saturates the inequalities (4-12) and turns the inequalities 
to equalities, that is, we have p- = -r^ , i = 2, 3, 4. The equation (4-15) indicates that Amax 

-'•■max ^ ' 

is a monotonic increasing function of p'^ and its maximum value corresponds to p'l = |, with 

A^„, = 2(l-pi). (4-16) 

and 

Substituting the results that obtained for A^^x ^-nd , i = 1, 2, 3, 4 in p — A^^xP's^ we obtain 

P - Amaxp's = (2pi - l)|0l)(0l|, (4-18) 

which is a pure state in agreement with theorem (2) of Ref.[3]. 
Therefore, we have 

p = 2(1 -pi)p1 + (2pi - l)\<Pi){<Pi\, (4-19) 
which is optimal LSD of BD states in agreement with [5] . 

4.2 Iso-concurrence decomposable states 

In this section we define iso-concurrence decomposable (ICD) states, then we obtain optimal 
LSD for this example. The iso-concurrence states are defined by [14, 15, 16] 

|0i) = cos^lOO) +sin^|ll)), 102) =sin^|00) -cos^lll)), (4-20) 
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103) = cos^lOl) +sin^|10)), 104) =sin^|01) -cos^llO)). (4-21) 

It is quite easy to see that the above states are orthogonal and thus span the Hilbert space of 
2 (g) 2 systems. Now we can define ICD states as 

p = ^p,|0,)(0,|, 0<p, <1, Y.Pi^^- (4-22) 

i=l i=l 

These states form a four simplex (tetrahedral) with its vertices defined by pi = 1, P2 = 1, 
Ps — 1 and P4 — 1, respectively. 

Peres-Horodeckis criterion [17, 18] for separabihty implies that the state given in Eq. (4-22) 
is separable if and only if the following inequalities are satisfied 



(Pl 


-P2) 


< 


^4p3P4/sin2 2^+ (p3 


-P4)^ 


(P2 


-Pi) 


< 


^4p3P4/sin2 2^ + (P3 


-P4)^ 


(P3 




< 


V4piP2/sin2 2e+{pi 


-P2)^ 


(P4 


-Ps) 


< 


^4pip2/sin2 2e + {pi 


-P2)^. 



(4-23) 
(4-24) 
(4-25) 
(4-26) 

Inequahties (4-23) to (4-26) divide tetrahedral of density matrices to five regions. The cen- 
tral regions, defined by the above inequalities, form a deformed octahedral and are separable 
states. In the other four regions one of the above inequality will not hold, therefore they 
represent entangled states. Bellow we consider entangled states corresponding to the violation 
of inequality (4-23) i.e. the states which satisfy the following inequality 



(Pi - P2) > V4P3P4/ sin^ 26 + (pa - Pa)^- (4-27) 
Now in order to obtain optimal LSD of entangled Iso-concurrence decomposable state given 



in (4-1), with (pi — P2) > \J^PzPa/ sin^ 29 + (p^ — Pa)^, we first choose an arbitrary separable 
state 



p's^J2p'i\<l>i) {<l>i\, 0<pi<l, X^P- = 1, p'l < P2 + V 4P3P1/ sin^ 20 + (p3 - p;)2. 



i=l 



(4-28) 
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in the separable region. Then using strict SDP optimization prescription of section (2), we try 
to optimize Tr{Ap'g) with respect to p — Ap^ > 0, where the feasible solution corresponds to 

^max = min{^, ^, ^, (4-29) 

Pi P2 P3 Pi 



Now, using the inequalities 



Amax<^, fori ^2,3, A (4-30) 
Pi 



and summing over the indices i=2,3 and 4, we obtain 

(l-p;)A„„, <(1-Pi), (4-31) 

since, we have 

^ > \^ ,Pi>P2 + V4p3P4/sin2 2^ + (p3-P4)2 > p[. (4-32) 
Pi 1 - Pi 

The only possible choice of Amax which is consistent with the positivity of p — AmaxP's is 

_ (1-Pi) . . 

(1 -Pi) 

This choice of A^ax given in (4-33) saturates the inequalities (4-30) and turns the inequalities 
to equalities, that is, we have p- = , i = 2, 3, 4. The equation (4-33) indicates that A^ax 

^ -'•■max ^ ' 

is a monotonic increasing function of p[ and its maximum value corresponds to 



P'l =P2 + ^^PzPj sin^ + (p3 - p;)2, (4-34) 

with 

Araax = \^ ■ (4-35) 
^ ~ Pi 



and 



= 7(^-^,^ = 2, 3, 4. (4-36) 
(1 - Pi) 



According to relation (4-34) we have 



Amax = 1 - (Pi - P2) + V4p3P4/ siu^ 29 + (ps " Pi? ■ (4-37) 
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Substituting the results that obtained for A^ax and , i = 1, 2, 3, 4 in p — Amaxp's-i we obtain 

/ C 

P - Amaxp's = ^^^l*^!)^*^!!' ^^^^^ 

( C is concurrence defined in [15]) which is a pure state in agreement with theorem (2) of 
Ref.[3]. 

Therefore, we have 



p = (1 - (pi -P2) + V4p3P4/sin2 2^ + (p3-P4)2)pl + (2pi - l)|0i)(0i|, (4-39) 

which is optimal LSD of ICD states in agreement with [5] . 

In the special case of {9 — 7r/4) we obtain Bell decomposable state and A^ax — 2(1 — pi). 

4.3 A generic 2x2 density matrix in Wootters's basis 

In this subsection we obtain optimal LSD for a generic two qubit density matrix by using 
Wootters basis. Wootters in [19] has shown that for any two qubit density matrix p there 
always exist a decomposition 

i 

called Wootters's basis, such that 

{xi\xj) = XiSij (4-41) 

where Aj are square roots of eigenvalues, in decreasing order, of the non-Hermitian matrix pp 
and 

p={ay(^ay)p*{ay(^ay) (4-42) 

where p* is the complex conjugate of p when it is expressed in a standard basis such as 
{| TT); I Ti)) I IT)) I ii)} and ay represent Pauli matrix in local basis {| j), | |)}. Based on 
this, the concurrence of the mixed state p is defined by max(0, Ai — A2 — A3 — A4) [19]. 
Now let us define states \x'j) as 

= /or i = 1,2, 3, 4. (4-43) 
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Then p can be expanded as 

p^j:Xi\x[){x[\ (4-44) 

i 

and Eq. (4-41) takes the following form 

WK) = % (4-45) 

Also Wootters has shown that p is separable if Ai — A2 — A3 — A4 < and if p is in boundary 
separable state then Ai = A2 + A3 -I- A4. By defining Pj = AjiCi, where ki — {x[\x'j), then 
normalization condition of p leads to 

Tr(p) = X:A = l. (4-46) 

i=l 

Now in order to obtain optimal LSD of entangled Wootters state given in (4-40), (4-41), with 
Ai > A2 + A3 -|- A4, we first choose an arbitrary separable state 

Ps = T.KWi){<\, o<X<r' T^Kki^i, a; < A'2 + A'3 + Ai. (4-4?) 

i=l i=\ 

in the separable region. Then using strict SDP optimization prescription of section (2), we try 
to optimize Tr(Ap'^) with respect to p — Ap'^ > 0, where the feasible solution corresponds to 

J^max = min{^, ^, ^, ^} = min{^, ^, ^, ^} (4-48) 
Ai A2 A3 A4 Pi P2 P3 Pi 



Now, using the inequalities 



A^„, <^, /ori = 2,3,4 (4-49) 
Pi 



and summing over the indices i=2,3 and 4, we obtain 

{l-p^)Kmax<{l-Pl) (4-50) 

since, we have 

^ > , Ai > A'2 + A'3 + A;> a;. (4-51) 
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The only possible choice of Amax which is consistent with the positivity of p — AmaxP's is 

''-"^-(i-p;)-(i-m;)- ^ ^ 

This choice of A^ax given in (4-52) saturates the inequalities (4-49) and turns the inequalities 
to equahties, that is, we have A' = , i = 2, 3, 4. The equation (4-52) indicates that Ajnax 
is a monotonic increasing function of A'^ and its maximum value corresponds to 

a; = A'2 + A'3 + A^, (4-53) 

with 

A - 

1 - X[ki ' 



1 — Ai fci , , 

Amax = T^, (4-54) 



and 



1 - Xih 

Therefore, we can show that 



a; = ^V\'T^'^ - (4-56) 



Amax = 1 - A:i(Ai - A2 - A3 - A4) = 1 - hC, (4-57) 

where C is concurrence. Using ki — {x\\x'j) and (4-43) we obtain 

C 

Amax^i- ^{Xl\xi). (4-58) 

M 

Substituting the results that obtained for Amax a-nd A^ , i = 1, 2, 3, 4 in p — AmaxP'sj we obtain 

P - AmaxP's = C\x[) {x[\, (4-59) 

which is a pure states in agreement with theorem (2) of Ref.[3] . 
Therefore, we have 

p = (1 - ^{x,\x,))p, + C\x[){x[\, (4-60) 
which is optimal LSD of states in agreement with [5] . 
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4.4 2 (g) 3 Bell decomposable state 

In this subsection we obtain optimal LSD for the Bell decomposable states of 2 (g) 3 quantum 
systems. A Bell decomposable density matrix acting on 2 (g) 3 Hilbert space can be defined by 

6 6 

P = Y.Pi\'^i){'^i\^ 0<pi<l, ^A = l, (4-61) 

i=i 1=1 

where \ipi) are Bell states in = ® Hilbert space, defined by: 





= i=(|ll) + |22)), 






- |22», 




= i=(|12) + |23)), 






- |23)), 




= -^(|13) + |21)), 




-7=2^^''^ 


-|21)). 



(4-62) 



It is quite easy to see that the above states are orthogonal and hence they can span the Hilbert 
space of 2(8)3 systems. Prom Peres- Horodeckis [17, 18] criterion for separability we deduce that 
the state given in Eq. (4-61) is separable if and only if the following inequahties are satisfied 

(pi -P2)' < (P3 +P4)(P5 +P6), (4-63) 

(P3 - PaY < (P5 + Pe) (pi + P2), (4-64) 

(P5 -Pef < {Pi + P2) (P3 + Pd- (4-65) 

In the following we always assume without loss of generality that Pi > P2, Ps > Pa and ps >Pq. 

Again in order to obtain optimal LSD of entangled BD state given in (4-61), with pi > 
P2 + (pa -l- Pa) (pb + Pe) , we first choose an arbitrary separable state 

p1 = EpI I'/'i) , 0<p^<l, E^i^l' p'i<P2 + V(^^3+P4)(P5 + P6);, (4-66) 

i=l i=l 

in the separable region. Then using strict SDP optimization prescription of section (2), we try 
to optimize Tr{Ap'g) with respect to p — Ap'^ > 0, where the feasible solution corresponds to 

» . rPl P2 P3 P4 P5 P6. «„x 

Amax^rnin{ — ,—,—}. (4-67) 
Pi P2 Pa P4 P5 Pe 
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Now, using the inequalities 

Ama. fori ^2,. ..,6 (4-68) 

Pi 

and summing over the indices i=2,... and 6, we obtain 

{l-p[)Amax<{l-pi), (4-69) 

since, we have 

^ > ^ ,Pl >P2 + ^iP3+PA){p5+P6) ,p'l < P2 + \/{P3+P4){P5+P6)- (4-70) 
Pi ^ ~ Pi 

The only possible choice of Amax which is consistent with the positivity of p — AmaxP's is 

A = (^-P^) (4-71) 

'-max /-I ' \ ■ \ / 

(1 - Pi) 

This choice of A^ax given in (4-71) saturates the inequalities (4-68) and turns the inequalities 
to equalities, that is, we have p[ = ,i = 2,..., 6. The equation (4-71) indicates that 

A^ax is a monotonic increasing function of p'l and its maximum value corresponds to p'l — 



P2 + V {P's + p'a) {P'b + p'g): with 



and 



Amax = =l-Pi+P2 + sJ{p^+Pa){P5+P&)- (4-72) 



p: = -^,z = 2,...,6. (4-73) 



Substituting the results that obtained for A^nax and , i = 1, 6 in p — AmaxP's-i we obtain 

P - Amaxp's = (1 - Amax)\(l>l){(t>l\, (4-74) 

which is a pure states in agreement with theorem (2) of Ref.[3]. 
Therefore, we have 



p = (1 -pi +P2 + \/{p3+P4){P5+P6))Ps + {Pl-P2 - V(P3 +P4)(P5 +P6))|0l)(0l|, (4-75) 

which is optimal LSD of 2 x 3, BD states. 
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The above choice of A^ax do not cover the whole set of separable states lying at boundary 
Pi — P2 + \J (Ps + Pi) (P5 + Pq) ■ Hence we should try other possible values of h.max ^ follows: 

PiKnax <Pi , i = 1, 4, 5, 6. (4-76) 

Summing over indices i = 1,2,..., 4 we obtain 

A < ^-P^-P^ u_77) 
1 - P2 - Pa 

Therefore, maximum feasible choice of A^ax is 

. _ 1-P2- P3 , . 

^^max -.111 V / 

which is possible if we choose 

Pi = ir^ 1,4,5,6, (4-79) 

■''-max 

Now substituting p'^ — ,i — 1, 4, 5, 6 in normahzation condition Tr{p'g) — Yh=i — ^ 

the separability equation (4-63) we can solve A^ax as a function of single variable p'^ and after 
optimizing it with respect to pg we get 

Amax = 1 - (P2 - Pi) - (P3 + P4) - ^(P5 + Pe), (4-80) 



and 



2pi-p5-p6 ' P5+P6-4p4 

P2 = ^X 'P'^^A ■ ^^"^^^ 



Similarly, following the above procedure for other possible choices, such as 



Amax = ^ ,i = 1,3,4,6, (4-82) 
Pi 



yields 



2pi - P3 - P4 ' P3 +P4 - 4p6 I. 
^^ = ^A '^^^^A ■ ^^^^^ 

'-max '-max 
Amax = 1 - (P2 - Pi) - (P5 + Pg) " |(P3 + P4)- (4-84) 

Of course, using this procedure we can obtain an optimal separate decomposition with rank-3 
entangled part of some particular given density matrices. 
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4.5 Werner states 

The Werner states are the only states that are invariant under local unitary operations. For 
d<^ d systems the Werner states are defined by [20] 

Pf = (('^ - + ('^^ - 1)^) ' -!</<!' (4-85) 

where / stands for identity operator and F — J2i,j |u) (j^l- ^ is shown that Werner state is 
separable iff < / < 1. 

Now in order to obtain optimal LSD of entangled Werner state given in (4-86), with —1 < 
/ < 0, we first choose an arbitrary separable state 

in the separable region. Then using strict SDP optimization prescription of section (2), we try 
to optimize Tr{Ap'^) with respect to p — Ap'^ > 0, where the feasible solution corresponds to 

A = m^M^l^ ii-^l = (4-87) 

The equation (4-93) indicates that Amax is a monotonic increasing function of /' and its 
maximum value corresponds to /' = 0, with 

Ama. = / + 1- (4-88) 

Substituting the results that obtained for A^ax in p — A^axP'sj we obtain 

F - I 

P - ^maxPs = P{f} - Amax(l + /)P{/=0} = /( ^2 _ (^-89) 

which is a pure states in agreement with theorem (2) of Ref.[3]. 
Therefore, we have 

pyy = (1 + f)p{f=o} + /( J^)- (4-90) 
which is optimal LSD of Werner states. 
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4.6 Isotropic states 

The d<S)d bipartite isotropic states are the only ones that are invariant under U®U* operations, 
where * denotes complex conjugation. The isotropic states old®d systems are defined by [21] 

PF = ]^ (/ - (V'+l) + F |V^+) (V^+l , < F < 1, (4-91) 

where = :^J2i is maximally entangled state. It is shown that isotropic state is 

separable when < F < | [21]. 

Now in order to obtain optimal LSD of entangled Isotropic state given in (4-92), with 
I < / < 1, we first choose an arbitrary separable state 

p's = Pf = (/ - IV'^) (V'+l) + F' |V^+) (V^+1 , < F' < 1, (4-92) 

in the separable region. Then using strict SDP optimization prescription of section (2), we try 
to optimize Tr{Ap'g) with respect to p — Ap'^ > 0, where the feasible solution corresponds to 

A_ = mi„{i|l,|i^} = ii^. (4-93) 

The equation (4-93) indicates that A^nax is a monotonic increasing function of F' and its 
maximum value corresponds to F' = ^, with 

A = - (4-94) 

■''-max -j^j • V / 

Substituting the results that obtained for A^ax in p — A^axP'sj we obtain 

Pf - AmaxPl/d = (1 - ^maxM'^){'ip^\ (4-95) 

which is a pure states in agreement with theorem (2) of Ref.[3]. 
Therefore, we have 

PF = ^^jrfl/d + (1 - Ama.)|^+)(^+| (4-96) 

which is optimal LSD of isotropic states. 
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4.7 One parameter 3 3 state 

Let us consider a one parameter state acting on = (g) Hilbert space as [13] 



2 



a 5 — a 



+ y(J+ + -^'J-, 2 < a < 5, (4-97) 
where 

IV'^) = 73(|11) + |22) + |33)), 

(7+ = I (|12) (12| |23) (23| + |31) (31|) , (4-98) 

a_ = |(|21)(21||32) (32| + |13)(13|). 
Pa is separable iff 2 < a < 3, it is bound entangled iff 3 < o; < 4 and it is distillable entangled 

state iff 4 < a < 5 [13]. 

Now in order to obtain optimal LSD of entangled one parameter 3(8)3 state given in (4-97), 

with 3 < a < 5, we first choose an arbitrary separable state 



2 

Ps = Pa' = ^ 



+ yf^+ + 2 < a' < 3, (4-99) 

in the separable region. Then using strict SDP optimization prescription of section (2), we try 

to optimize Tr{kp'g) with respect to p — Ap'^ > 0, where the feasible solution corresponds to 

5 — oi 

Amax = 7 r (4-100) 

The equation (4-100) indicates that Aj^ax is a monotonic increasing function of a' and its 
maximum value corresponds to a' = 3, with 

Amax = (4-101) 

Substituting the results that obtained for A^ax in p — AmaxP's^ we obtain 

P - AmaxPa=3 = (2/7|V'+)(V'+| + 5/7(7+) (4-102) 

which is a pure states in agreement with theorem (2) of Ref.[3]. 
Therefore, we have 

5 — Oi 

P = {^)Pa=3 + (2/7|V^+)(V^+| + 5/7a+) (4-103) 
which is optimal LSD of one parameter 3 3 states. 
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4.8 Multi partite isotropic states 

In this subsection we obtain the optimal LSD for a n-partite d-levels system. Let us consider 
the following mixture of the completely random state po — I/d^ and the maximally entangled 
state IV'"'") 

p{s) = (1 - + 5 U+ , < s < 1, (4-104) 



where / denotes identity operator in ^"-dimensional Hilbert space and j?/'"'") = Sf=i \ii • • ■ i)- 
The separability properties of the state (4-104) is considered in Ref. [22]. It is shown that the 
above state is separable iff s < sq = (1 + o?"~^)~^. 

Now in order to obtain optimal LSD of entangled multi partite isotropic state given in 
(4-104), with So < s < 1, we first choose an arbitrary separable state 

= p{s') = (1 - 5')4 + i^^) (i^^ , 0<s'< So, (4-105) 



0?" 

in the separable region. Then using strict SDP optimization prescription of section (2), we try 
to optimize Tr(Ap'^) with respect to p — Ap'^ > 0, where the feasible solution corresponds to 

Ama. = (4-106) 

I — s' 

The equation (4-106) indicates that A^ax is a monotonic increasing function of s' and its 
maximum value corresponds to s' — Sq, with 

Substituting the results that obtained for A^ax in p — Aj^axP's-i obtain 

P{S) - AmaxP{s = So) = (1 - AmaxW){i^^\ (4-108) 

which is a pure states in agreement with theorem (2) of Ref. [3]. 
Therefore, we have 

which is optimal LSD of multi partite isotropic states. 
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Here in this work we have been able to obtain LSD of bunch of mixed state density matrices via 
an elegant method of convex positive semidefinite optimization methods, where the results that 
obtained are in agreement with those obtained by other methods in Ref. [14, 15]. Comparing 
this method with those of previously introduced one, one can appreciate the elegance and 
usefulness of SDP method in connection with LSD. 
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